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ABSTRACT 

We present a catalogue of 348 galaxy clusters and groups with 0.2 < z < 1.2 
selected in the 2.78 deg^ ALHAMBRA Survey. The high precision of our photometric 
redshifts, close to 1%, and the wide spread of the seven ALHAMBRA pointings ensure 
that this catalogue has better mass sensitivity and is less affected by cosmic variance 
than comparable samples. 

The detection has been carried out with the Bayesian Cluster Finder (BCF), whose 
performance has been checked in ALHAMBRA-like light-cone mock catalogues. Great 
care has been taken to ensure that the observable properties of the mocks photometry 
accurately correspond to those of real catalogues. From our simulations, we expect 
to detect galaxy clusters and groups with both 70% completeness and purity down 
to dark matter halo masses of Mh 3 x IO^^Mq for z < 0.85. Cluster redshifts are 
expected to be recovered with 0.6% precision for z < 1. We also expect to measure 
cluster masses with 0.25 — 0.35dea; precision down to ^ 3 x IO^^Mq, 

masses which are 50% smaller than those reached by similar work. 

We have compared these detections with previous optical, spectroscopic and X- 
rays work, finding an excellent agreement with the rates reported from the simulations. 
We have also explored the overall properties of these detections such as the presence 
of a colour-magnitude relation, the evolution of the photometric blue fraction and 
the clustering of these sources in the different ALHAMBRA fields. Despite the small 
numbers, we observe tentative evidence that, for a fixed stellar mass, the environment 
is playing a crucial role at lower redshifts (z<0.5). 

Key words: catalogues, cosmology: large-scale structure of Universe, cosmology: ob¬ 
servations, galaxies: clusters: general, galaxies: clusters: individual, galaxies: evolution 


1 INTRODUCTION 

Galaxy clusters are the largest objects gravitationally bound 
in the universe. According to the standard model of cosmic 
structure formation, they appeared as a result of the initial 
perturbations in the mass power spectrum at a typical co- 

* E-mail: begona.ascaso@obspm.fr 


moving scale of ~ 10h~^ Mpc. At larger scales, the universe 
is mainly dominated by gravity with the gas dynamics a mi¬ 
nor contributor. However, at smaller scales, the complexity 
of the astrophysical processes, particularly related with for¬ 
mation and evolution of galaxies, produces changes in the 
observational properties of the structures. 

The first catalog ues of galaxy clusters came by th e hand 
of Abell (lAbelllll95a ) and Zwicky (IZwickv et al.l[l96l]) in the 
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early sixties, together with po sterior improved e xtensions ex¬ 
panded to the southern sky (I Abell et al.lll98^ . These cat- 
alognes, in addition to snffering from large projection ef¬ 
fects due to the absence of multi-band photometry or spec¬ 
troscopy, had a complicated selec tion fnnction since they 
were based on visual inspections JKatgert et al.lfl99fil and 
references herein). 

Snbsequent cluster and group catalogues have been 
bnilt from systematic searches of galaxy clusters in 
different wavelengths or with different techniques: (1) 
searches in optical data (for a revi ew of th e differ¬ 
ent m ethodologies in the optical, see iGall [200^ Ascas^ 
201S); (2) searches with X-ray data, (see iRosati et al.l 


2002 ; iBurenin et al.l l2007l and references herein) ; (3) 
searches nsing radio sources fe.g. ICalametz et al.l 120091: 
Chiaberge et al.|l2010l:[Castignani et al.l 20141: iBlanton et al.l 


201A) (A) searches using the Sunyaev-Zel’dovich (SZ) sig- 


natnre in cosmic microwave b ackground maps (e.g. iBartletu 
l2004l : [^aniszewski et al.ll2009h: (5) searches using the Weak- 
Lensing (WL) effect (see Wittman et al.l l2006l for a de¬ 
tailed explanation and ref erences herein); (6) spect roscopic 
searches (see, for instance iKnobel et al.ll2009l . I20l3 and ref¬ 
erences arising from them). All these different techniques 
provide well-characterized selection functions, completeness 
and purity rates. 

Since the beginning of the XXI century, we have been 
witnesses of the discovery of even more extreme clus¬ 
ters. For the first time, high redshift {z > 1) struc¬ 
tures in the unive rse using optical/IR data, X-rays, radio 
or th e SZ effect |lGalametz_et_^ [200^; iFassbender et aP 
201ll: IPlanck Collaboration et J] l201ll: Jee et all l201ll : 


Brodwin et al. I I2OI3I: ICastignani et al.ll2014l . to name a few) 
have been found. The finding of these clusters, usually very 
massive (several times IO'^'^Mq), has already impulsed a 
change in the main theories of cluster formation and galaxy 
evolution. For instance, the discovery of very ma ssive clus- 
ters such as XMMU J2235.3 - 2 557 at z ^1.4 (I Jee et al.l 
l2009h or ‘El Gordo’ at ~0.87 (iMenanteau et ^ 201^ 
has c hallenged the main cosmological theories i Hovle et al.l 
l2012h . The hrst proto-clusters at redshift > 2 have also 
set u pper limits in the time scale of struc ture formation 
(e.g. ICapak et al.l I2OIII : IChiang et al.l l2014h . Complemen- 
tarily, the di scovery of post-me rger clusters (e.g. the ‘Bul- 
let Cluster’, ^adac et al.ll2008l . the ‘Musket Ball cluster’, 
iDawson et al. ll2Q]J) has shown that clusters are far from 
being static entities, being able to merge and dramatically 
change their properties. 

The observational strategies followed to construct the 
present cluster surveys have favored the detection of the 
most massive and luminous clusters due to fl ux limits and 
resolution e ffects. Only spectroscopic searchesjKnobel_et_^ 


200S, I 2 OI 2 II and recent X-ray surveys (iFinoguenov et ^ 


20071), have provided complete group catalogues down to 
low masses (~ IO'^^Mq). Unfortunately, the observational 
cost of these surveys is very high and only small areas have 
been sampled. 

In the last few decades, several multiple medium- 
bands surveys have been developed: the Classifying Ob¬ 
jects by Med ium-Band Obser vations in 17 Filters survey 
fGOMBO-17. IWolf et al.l 2003); the Cosmic Evolution Sur¬ 
vey fCOSMOS. IScoville et al.l 2007al '): the Advanced Large, 
Homogeneous Area Medium Band Redshift Astronomical 


survey (ALHAMBRA, iMoles et al.l l2008lj : or the Survey 
for High-z Absorption Red and Dead Sources (SHARDS, 
iPerez-Gonzalez et al.l[^13lj . among others. One of the main 
benefits of having a large number of medium-band filters 
covering the whole optical spectrum at least is the fact that 
the photometric redshift resolution becomes comparab le to 
that of a spectroscopic survey fsee lMolino et aHl2014l for a 
review) allowing us to sample low-massive structures in the 
universe within larger areas than spectroscopic samples. 

In this work, we have used the ALHAMBRA survey to 
perform a systematic search of galaxy cluster s and groups 
using the Bayesian Gluster Finder (BCF, lAscaso et al.l 
l20ial . The ALHAMBRA survey consists of four square de¬ 
grees imaged in 20 optical narrow bands and three broad¬ 
band IR bands. For technical reasons, only three degrees 
were observed and calibrated. In addition, part of the images 
were masked to take into account image artifacts and satu¬ 
rated stars, resulting into a final usable area of 2.78 square 
degrees. The survey is complete down to F’8 14IU < 24.5, 
wher e 7^81414^ is a synthetic combined band dMolino et al.l 
l2014h equivalent to the band with the same name at the 
Hubble Space Telescope (HST). The overall photometric 
redshi ft accuracy obtaine d for the survey is Az/)! -|- Zs) < 
0.015 dMolino et al.l[2oT3j . making it comparable to low res¬ 
olution spectra for each object in the survey. This survey, to¬ 
gether with its future ‘big brother’, the Javalambre- Physics 
of the accelerated universe Astrophysical Survey (J-PAS, 
iBenitez et al.l 1^141 1 will be able to set a benchmark in the 
determination of the cluster mass function in surveys. 

The structure of the paper is as follows. In section §2, 
we describe the ALHAMBRA dataset. Section §3 provides 
the basic information about the BGF cluster finder. Sec¬ 
tion §4 is devoted to the description of the mock catalogues 
used in this work on one hand and, on the other, to the 
results regarding cluster detection on them, the description 
of the selection function for the ALHAMBRA survey and 
the accuracy in measuring cluster properties such as red¬ 
shift or mass. In Section §5, we present the ALHAMBRA 
optical cluster and group detections. We first compare the 
detections found in this work with those found by other 
authors using different datasets and techniques. Then, we 
explore the main properties of the detections such as the 
stellar mass distribution, the presence of colour-magnitude 
relations and the fraction of blue and late-type galaxies in 
the cluster. Finally, section §6 provides the main conclusions 
of the paper and includes a discussion of the main results of 
this work. Where appropriate, we have used Ho=73 km s”*^ 
Mpc”*^, Hm =0.25, Ha=0.75 throughout this paper in or¬ 
der to match the same cosmology than the mock catalogues 
utilized in this work. 


2 THE ALHAMBRA SURVEY 

The Advanced Large, Homogeneous A rea Medium Band 
Redshift Astronomical (ALHAMBRAQ, iMoles et al.l l2008lj 
survey is a 20 narrow-band optical and three broad-band 
NIR (JHK) photometric survey imaged with the wide- 


^ http://alhambrasurvey.com/ 
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field cameras LAICA and OMEGA-2000 respectively, at the 
Calar Alto Observatory (Spainj^- 

This survey covers 4 square degrees spread in eight dif¬ 
ferent regions of the sky down to ~ 24.5 AB magnitude 
in the synthetic combined FSl^W band and ~ 20.5 in the 
infrared bands. Due to technical issues related with the un¬ 
availability of calibration stars in one of the fields (ALH-1) 
and the lack of two pointings in the ALH-4 and ALH-5 fields, 
only three degrees were available. Also, additional areas of 
the images such as saturated stars and image edges with 
insufficient exposure time to provide accurate photometry 
have been prop erly masked, obtai ning a final area of 2.78 
square degrees (iMolino et al.ll2014l ~l. 

The different fields of the ALHAMBRA survey were 
chosen strategically in order to overlap with well-known 
fields, many of them with multiwavelength data avail¬ 
able such as the Sloan Digi tal Sky Survey (SDSS, 


lYork et al.ll200d : lAhn et al.ll2014 L the DEEP2 Galaxy Red- 
shift Survey (iNewman et al.l 2013 1 or the GOSMOS survey 
dScoville et ah 2007a 1. 

The photometric redshifts obtained for this survey have 
been calculated with BPZ2.0 (Benitez in prep), an improved 
version of BPZ llBemtej|2000h capable of also providing stel¬ 
lar masses. The catalogues, including a full range of mea- 
sured properties of the galaxies, are publicly availabl^l (see 
iMolino et al .1120141 for details). Additionally, a set of mask 
files including saturated star and spurious effects are also 
available for the survey. These were generated together with 
the moc k catalogues and t heir building procedure can be 
found in IMolino et al.l[2014l . 


3 THE BAYESIAN CLUSTER FINDER 

In this work, we have us ed the Bayesian Gluster Finder 
(BGF; lAscaso et al. lEoIi), a technique developed to de¬ 
tect galaxy clusters and groups in any optical/infrared im¬ 
age dataset. This method i s based on a modifica tion of the 
Matched Filter Technique jPostman et al.ll2002li including 
photometric redshifts and the presence of Bayesian priors. 

In more detail, the BCF initially calculates the prob¬ 
ability that there is a cluster centered on each galaxy at a 
given redshift. In order to calculate this likelihood, we as¬ 
sume that clusters are modeled as a convolution of a partic¬ 
ular density, luminosity and photometric red shift profile. We 
choos e to use the Plummer density profile (IPostman et al.l 
l2002h . the Schechter luminosity function ([Schecht^ 1976) 
and the full redshift p robability function fPfzL iBenited 
I 2 OOOI: IMolino et al.ll20l4 in case it is available or a Gaussian 
approximation otherwise. This likelihood probability does 
not include any pre-assumption about the colours of the 
cluster. We convolve this likelihood with a prior probabil¬ 
ity to obtain the final probability. The prior refers to those 
properties that are not necessarily present in all clusters but 
can help to discern between different solutions. We choose 
to model two main features that are present in the major¬ 
ity of the clusters up to redshift ~ 1.6 and down to masses 
lO^^M© at least: the presence of the colour-magnitude re- 
lation (CMR; e.g. iLopez-Cruz et al.l [20041 : iMei et aklliood : 


lAscaso et al.l I 2 OO 8 I : IPapovich et al.l l201Cll~l and the presence 
of a well-defined brightest cluster galaxy (BCG), following a 
tight relati on between its magnitude a nd the redshift of the 
cluster (e.g lAscaso et al.ll201ll . l2014bll . 

In order to characterize the former prior component, 
we first created synthetic g — i and i — z colours predic- 
tions using a typical elliptical spectru m from the lib r ary b y 
IColeman et al.l lll980l'l as performed in lAscaso et al.l (I2OI2II . 
We chose those bands in order to be able to sample the 
4 OOOA break at low (2 < 0.7) and high (2 > 0.7)-redshift 
range respectively. We artificially created these bands by 
calculating the contribution of each of the ALHAMBRA 
narrow band s to th e new synthetic band as performed in 
IMolino et al.l (l2014l'l . We then created a Gaussian filter char¬ 
acterizing the colour for each redshift slice. The width of this 
Gaussian has been set to 0.5 in order to account for larger 
dispersions in the RS. As for the latter prior feature, we 
first measured empirically the K-band Hubbl e diagram for 
a com p lete sample of BC Gs e xtracted from lYTijev et al.l 
1I2OO8II . IStott et ^ 1I2OO8II and ICollins et al.l 1 2009l ~) up to 
redshift < 1. Then, we performed a colour transformation 
to our reference band, ESldlF using the same library of syn¬ 
thetic templates as for the CMR. Similar to the case of the 
CMR, we created a Gaussian filter characterizing the mag¬ 
nitude of the BCG at each redshift slice. Note that for a real 
detection, the amplitude of the likelihood is always several 
orders of magnitude larger than the signal of the prior. As a 
consequence, the prior information helps to discern between 
different solutions at different redshift slices without losing 
preliminary detections by the original likelihood. In other 
words, the BCF does not specifically rely on the colours of 
the galaxies to select clusters, therefore it is able to detect 
any structure over the threshold limit, independent of its 
colour. 

We performed a search in a predefined number of red¬ 
shift slices. The minimum redshift threshold (zmin) comes 
from the angular extent of the survey which is limited by its 
geometry; and the maximum redshift (zmax) is estimated 
from the wavelength coverage and the depth of the survey. 
The bin width (zbin) is fixed to be three times the expected 
photometric resolution of the survey. For instance, for the 
ALHAMBRA survey, we have fixed Zmin = 0.2, Zmax = 1.2 
and Zbin = 0.05. Effects of star s masking and edge s of the 
frames are taken into acco unt (IMolino et al.l l2014l l . As in 
lAscaso et al.l (|2012| . l2014al L we applied a probability cor¬ 
rection proportional to the effective area within 0.5 Mpc of 
the considered galaxy. This correction accounts for galaxies 
lying close to the border of the image or saturated stars. 

Afterward, clusters are selected as the density peaks 
of those probability maps and the center is located at the 
peak of the probability. Finally, if we find a cluster or group 
detected in different redshift slices (two or more detections 
separated by less than 0.5 Mpc and one redshift bin dif¬ 
ference), we merge them into one. For a mor e detailed de¬ 
script i on of t he method, we refer the reader to lAscaso et al.l 
Il20i2l.l2014'i) . 

The output of the algorithm consists of the centered po¬ 
sition, a measurement of the redshift and a measurement of 
the richness of the cluster. We have used the Acl parameter, 
defined as 


^ http://www.caha.es/ 

® https://cloud.iaa.csic.es/alhambra/ 


Acl = 


J2tiH{R<0.5Mpc) 


L* 


( 1 ) 
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i.e., the sum of the luminosity of the galaxies (in the 
-F814VK band) statistically belonging to the cluster divide d 
by the characteristic luminosity L* (jPostman et al.ll2002h . 
The galactic population of the cluster is defined as those 
galaxies lying within a given radius (0.5 Mpc, for the rich¬ 
ness calculation) with a given cut in photometric redshift 
odds {odds> 0.5/(l -I- z), in this work) and that the differ¬ 
ence with their photometric redshift and the redshift of the 
cluster is 


\z — Zch\ < 2bin/2 -I- (JnMAd(1 + 2Cl) 


( 2 ) 


where z is the redshift of the galaxy, zcl is the redshift of 
the cluster, Zhin is the redshift bin used for the detection 
and ctnmad is the expected photometric redshift accuracy 
of the survey (0.0125 for the case of ALHAMBRA). The 
latter parameter, (Jnmad, is defined as 


crNMAD = 1.48 X median 


\d- <d>\ \ 
1 + Zs J 


( 3 ) 


and d = zcl — Zs (see, for instance iMolino et al.l[20l3 L 

In this work, we have introduced a new richness mea¬ 
surement in addition to the Acl parameter, the cluster total 
stellar mass, Mql, 


N 

Mcl. = ^ M* (R < 0.5Mpc) (4) 

i=l 


i.e, defined as the sum of all the stellar masses of the galax¬ 
ies, M*, statistically belonging to the cluster. The stellar 
m asses have been calc ulated with BPZ2.0 in the same way 
as iMolino et al.l (l2014l L 

So far, the BCF has been applied to two more op¬ 
tical surveys: a wide s urvey, the CFHTL S -Archive Re- 
search Survey (CARS, lErben et al.l l200ill : lAscaso et al.l 
2012 h. and a very de e p survey, the Deep L ens Survey (DLS, 
Wittman et al. l2002l : lAscaso et al.ll2014^ L In this work, we 
apply the BCF to a high photometric redshift resolution 
survey, the ALHAMBRA survey. The comparison between 
these studies will indicate the benefits and drawbacks of us¬ 
ing datasets with different properties for detecting galaxy 
clusters and groups in the optical. 


4 SIMULATIONS 

In this section, we first describe the light-cone mock cata¬ 
logues that we have used to mimic the ALHAMBRA data 
(§4.1). Then, we use the BCF to detect galaxy clusters and 
groups in those mocks and obtain accurate cluster and group 
selection functions (§4.2) and we finally explore the dark- 
matter halo mass-richness relation in §4.3. 


4.1 Light-cone mock catalogue 

The light-cone mock catalogue that we have util i zed in 
this anal ysis has been obtained f rom iMerson et al.l (l2013l l 
(see also lArnalte-Mur et al.lf2014h . This mock galaxy cat¬ 
alogue has been built from a semi-analytical model of 
galaxy formation, applied to the halo merger trees extracted 
from a cosmologica l N-body simulatio n. The semi-analytical 
model used is the iLagos et al.l (1201 il l variant o f the semi- 
analytical galaxy formation model GALFORM (ICole et al.l 


l200(]ll . which models the star formation and merger history 
for a galaxy. Among other physical processes, this model 
includes feedback as a result of SNe, active galactic nuclei 
(AGN) and photo-ionization of the intergalactic medium. 
The model predicts the star formation history of the galaxy 
and therefore its spectral energy distribution (SED). The 
population of dark matter (DM) haloes for t he mock cata¬ 
logue comes from the Millennium Simulation (ISoringel et al.l 
12(30511 ■ a 2160^ particle N-body simulation of the A Cold 
Dark Matter cosmology starting at z = 127 and models the 
hierarchical growth to the present day. The halo merger trees 
are constructed using particle and halo data stored at 64 
fixed epoch snapshots spaced logarithmically in expansion 
factor. The minimum halo resolution is 20 particles, corre¬ 
sponding to 1.72 X 10^°/i“^Mq. Finally, the light-cone was 
constructed from this simulation by replicating the simula¬ 
tion box and choosing an orientation. In addition, a flux cut 
at FSIAW < 24.5 AB was applied to reproduce the selec¬ 
tion of the ALHAMBRA survey. The final mock catalogue 
is limited to z < 2 an d it does not includ e stars. All the 
details can be found in Merson_et_aljj 2013l3. 

As performed in lArnalte-Mur et al.l (I2014l3 . we created 
fifty non-overlapping realizations of the ALHAMBRA sur¬ 
vey mimicking its geometry from the whole mock catalogue. 
Besides, we manually included saturated stars and edge ef¬ 
fects in the mock catalogues resembling the masks extracted 
from the ALHAMBRA data. 

Initially, we ran the well-known photometric redshift 
code, the Bayesian Photometric Redshift (BPZ2.0, iBenfted 
120001 Benitez in prep). This code, the same used to calcu¬ 
late the ALHAMBRA photo-z, has a library of empirical 
templates with a very low outli er rate in high qua lity photo¬ 
metric catalogues (< 1% in the lllbert et al.lf2009l COSMOS 
catalogue, Benitez, private communication; < 2% in the AL¬ 
HAMBRA data, and most of those seem to be stars or active 
galactic nuclei (AGN), see IMolino et al.ll2014l3 . 

Since photometric redshifts are exquisitely sensitive to 
any discrepancy between the template library and actual 
galaxy photometry, it can be concluded that the rather 
sparse (11) BPZ2.0 library contains a complete (up to a 
few %), even if coarse-grained, representation of real galaxy 
colours for the galaxy populations sampled by the intersec¬ 
tion of ALHAMBRA and the spectroscopic redshift cata¬ 
logues used to measure the outlier rate. Therefore, any sys¬ 
tematic mismatches between galaxy colours in mocks and 
the BPZ2.0 template library signal the presence of non- 
realistic galaxy types in those mocks. Note than non-realistic 
does not mean that those types are physically absurd; only 
that its aggregated frequency nature is significantly below 
the 1 — 2% outlier rate we observe in the real world. 

This is exactly what we find when we run BPZ2.0 on 
the original mock catalogue, obtaining a photometric red¬ 
shift accuracy of Az/( 1 -|- z) = 0.0319, a f actor of ~ 3 worse 
than in the real data (IMolino et al.l[20l33 . A closer inspec¬ 
tion of the original photometry included in the light-cone, 
evidenced the above mentioned fact that the SED of a sig¬ 
nificant fraction of the galaxies included in this mock did 
not match any of the spectra in the BPZ.2.0 library. There¬ 
fore, to obtain a galaxy colour distribution which accurately 
resembles our real ALHAMBRA catalogues, we had to force 
these objects into realistic colours, assigning then the closest 
BPZ2.0 SED. This technique, called PhotReal, has already 





























































Galaxy clusters and groups in the ALHAMBRA Survey 5 



Figure 1. Density plots of the photometric redshift versus spec¬ 
troscopic redshift for the initial mock catalogue (left panel) and fi¬ 
nal photometry corrected mock catalogue (right panel). The den¬ 
sity scale is shown in the right part of the plot and the overall 
photometric redshift resolution (ctnmad) is quoted (see Eq. [Sjl. 
The photometry correction clearly improves the quality of the 
photometric redshift and make them comparable to the data. 

been applied to prev i ous work (lArnal te-Mur et al.l I20l4 
IZandivarez et af]l2014l : lAscaso et al.ll2015l ). will be fully de¬ 
tailed in a future publication (Benitez et al. in prep). 

Once the SED is thus chosen, we generated galaxy fluxes 
through the ALHAMBRA set of filters and add to them em¬ 
pirically calibrated photometric noise. This noise is a combi¬ 
nation of the expected photometric noise from the observed 
relationship between magnitudes and errors in the ALHAM¬ 
BRA filters, plus a systematic noise which is approximately 
constant with magnitude and most likely unavoidable when 
measuring galaxy colours in multiband photometry. This 
systematic is empirically calibrated to be 8% for bluer ob¬ 
jects and 6% for red galaxies (Benitez, in preparation). 

Afterwards we run BPZ2.0 on those mock catalogues. 
In Fig. m we display the density maps of the photometric 
redshift distribution versus spectroscopic redshift for the ini¬ 
tial mock catalogue (left panel) and the final mock catalogue 
(right panel) together with the overall photometric redshift 
resolution. We notice how, after performing this technique, 
we increase the resolution by a factor of ~ 3, reaching val¬ 
ues which are very similar to those of real data. The fact 
that the real and simulated photometric redshift values are 
so similar strongly support the fact that the BPZ2.0 library 
is a faithful representation of the real galaxies. 

This procedure of course implicitly introduces an error: 
there may be real galaxies which are present in the cat¬ 
alogues, but not covered by the BPZ2.0 templates or not 
adequately present in the spectroscopic redshift catalogues 
used to measure the outlier rate. Since we are not aware 
of any such substantial population within the ALHAMBRA 
depth, we estimate that the procedure followed here may 
introduce at most a contamination of a few percent, much 
smaller than that produced by o t her so urces. 

As stressed in iMolino et al.l ll20l4l with the ALHAM¬ 
BRA data, we can additionally increase this resolution by 
performing a cut in the odds parameter. This parameter is 
defined as the integral of the redshift probability function 
P{z) around its maximum peak within twice the expected 
photo-z accuracy for the survey/data which, in the case of 
the ALHAMBRA data, it was measured to be 0.0125. The 


odds parameter gives us a direct estimation of the quality 
of the photometric redshift since it measures how concen- 
tra ted around the ‘true’ photometric r e dshift value the P(z ) 
is llBenftezI 120001 : iBenitez et al.l l2009l : iMolino et al.l l20ld ) . 
For instance, performing a cut in odds > 0.5 in the simula¬ 
tion increases the photo-z resolution to l\z/{l-\-z) = 0.0098 
while a cut in odds >0.9 provides a subsample of even higher 
photometric resolution, A2/(l -|- 2) = 0.0047. Of course, a 
precise understanding of the selection function that the odds 
parameter introduces needs to be quantified. However, this 
parameter becomes very useful when detecting galaxy clus¬ 
ters for two main reasons: first, it decreases the size of the 
sample but still leaves ‘useful’ galaxies to detect galaxy clus¬ 
ters and second, we get rid of the field galaxy contamination 
which have usually lower odds parameter values on average. 
We will come back to this later in the section. 

Furthermore, we have performed two additional checks 
for this mock catalogue. The first one has been the compar¬ 
ison of the number counts per square degree between the 
mock catalogue and the real data as shown in Fig. (2] As 
we see here, the mock catalogue FSldVF synthetic magni¬ 
tude distribution traces almost perfectly the observed one 
once we removed the effect of stars by making a cut in the 
stellar flag < 0.7. The restriction of the ALHAMBRA data 
to 2 < 2 is plotted in order to match the redshift limit of 
the mock catalogue. The simulation counts are only ~ 8% 
smaller than the real data, a di screpancy that could be at - 
tributed to cosmic variance (e.g. lLopez-Saniuan et al']|2014l ). 

The second check performed has been to compare the 
stellar masses provided by BPZ2.0 in the mock catalogue 
with the stellar masses in the ALHAMBRA data. In Fig. 
O we display the stellar mass histograms for the ALHAM¬ 
BRA data, together with the initial simulation and the post- 
processed simulation. The stellar masses of the initial simu¬ 
lation refer to those provided by GALFORM. The three dis¬ 
tributions have been restricted to FSldlF < 24.5, the com¬ 
pleteness limit for the ALHAMBRA survey. We find that 
the agreement of the observed and post-processed stellar 
mass distribution is very good (<8%), while the difference 
is higher for the observed and initial stellar mass distribution 
(<17%). This l atter difference has already been reported in 
iMitchell et al.l ll2013l) when comparing the GALFORM stel¬ 
lar masses with stellar masses estimated using SED fitting, 
and could be due to the dust extinction applied to the model 
galaxies, which can be substantial for massive galaxies. The 
disagreement found at the low-mass end of the distribution 
can be justified as both, the survey and the mock catalogues 
are flux-limited. 

4.2 Detecting clusters in the mock catalogue 

We have considered as the reference sample in the simulation 
those haloes more massive than Mh ^ IO'^^Mq, together 
with the galaxies associated to each halo and the position of 
the center, set in one of the galaxies. Then, we have used the 
BGF, as described in section §3, to search for galaxy clusters 
and groups in the ALHAMBRA mock catalogue restricted 
to odds >0.5/(1-|-2). This odds cut has been performed in 
order to use the best photometric redshift quality galaxies 
and keep a constant galaxy density as a function of redshift. 

We have run the BGF on a test sample using both red¬ 
shift filters: the full redshift probability function (PDZ) and 
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Figure 2. F%liW magnitude distribution for the ALHAMBRA 
data and the mock catalogue data. The dotted-long-dashed black 
line and the green dotted-short-dashed line refers to the whole 
ALHAMBRA sample including and excluding the stars, respec¬ 
tively. The black solid line, refers to the ALHAMBRA sample 
excluding the stars up to redshift 2. The red dotted line and blue 
dashed line displays the F81AW magnitude distribution for the 
original and final mock catalogues, respectively. We find a very 
good agreement (<8%) between the FSliW distribution for the 
simulation and the ALHAMBRA data sample up to redshift 2, 
confirming the accuracy of the mock catalogue. 


a simple Gaussian centered in the redshift of the cluster. In 
this case, the results are very similar (<5% difference). Al¬ 
though a simple Gaussian is not the best approach in pho¬ 
tome tric surveys (iLopez-Saniuan et al.l 120141 : iMolino et al.l 
l2014h . the bright galaxies centered in clusters present high 
S/N ratios, making reasonable the approximation. There¬ 
fore, we decided to use the second approach for computa¬ 
tional and disk space purposes. 

We have measured the rate of completeness, defined as 
the percentage of clusters detected out of the total simulated 
sample, and the purity, defined as the percentage of clusters 
simulated that were detected out of the total detected sam¬ 
ple. In the bottom and top panel of Fig. [4] we show these 
rates as a function of dark matter halo {Mh) and total stel¬ 
lar mass (Mcl)) respectively. The dark matter-total stellar 
mass relation has been directly calibrated from the simula¬ 
tion and we investigate it in the next subsection. 

While the purity rates remain >60% for all redshift and 
mass richness, the completeness rates are lower and decrease 
with redshift in general. We obtain completeness and pu¬ 
rity rates >80% for clusters with total stellar masses larger 
than 1.5 x lO^^M© or halo masses of 6 x lO^^M© up to red¬ 
shift 0.8. For higher redshifts (0.8 < z ^ 1), we increase 
the threshold mass for which we find completeness >80% 
to Mql > 1-8 X IO^^Mq, equivalently Mh > 1.0 x lO^^M©. 
We fail at detecting galaxy clusters with completeness and 
purity rates >80% for clusters at redshift z > 1. 

If we relax instead both completeness and purity rates 



Figure 3. Stellar mass distribution for the ALHAMBRA data 
and the mock catalogue. The black solid line refers to the whole 
ALHAMBRA excluding the stars up to redshift 2. The red dot¬ 
ted line alludes to the GALFORM stellar mass distribution in 
the original mock catalogue and the blue dashed line traces the 
stellar mass distribution for the post-processed mock catalogue. 
All the distributions are restricted to the magnitude limit where 
ALHAMBRA is complete (F814VF < 24.5). We confirm a good 
agreement between the stellar mass distribution of the ALHAM¬ 
BRA data and the final mock catalogue (<8%), finding a higher 
difference (< 17%) between the ALHAMBRA data and the initial 
mock catalogue. 


to be higher than 70%, we obtain a stellar mass limit of 
1.2 X IO'^^Mq or halo masses of 3 x up to redshift ~ 

0.85, increasing to 1.5x lO^^M© or halo masses of 6x lO^^M© 
between 0.85 < z ^ 0.95 and to 1.8 x lO^^M©, equivalently 
1.0 X 10^'‘Mq, within 0.95 < 2 ^ 1.2. 

We have examined the same rates as a function of red¬ 
shift for different Acl ranges, as shown in Fig.[3 As before, 
the Acl to dark matter calibration has also been measured 
from the mock catalogue. We see a very similar behaviour 
of the purity rates, being this >60% for all redshift and 
mass ranges. Both completeness and purity become >80% 
for Acl > 51.5 up to redshift 0.8 and for Acl > 60.8 for the 
redshift range 0.8 <2^1. 

Decreasing the completeness and purity rates to 70%, 
results into a lower Acl threshold. We will be able to detect 
galaxy groups down to Acl > 41 up to redshift 0.8, down 
to Acl > 51.5 between redshift 0.85 and 0.95 and down to 
Acl > 60.8 for redshifts higher than 1. 

We have then summarized the selection function ob¬ 
tained from this analysis with three levels collected in Table 
[T] The first level, includes the mass and redshift limits for 
our detections to be at least 80% complete and pure and the 
second level refers to a level of, at least, 70% of completeness 
and purity. We have also included a third level, of those de¬ 
tections that have a level of completeness higher than 50% 
and purity higher than 60%. 

We have also investigated the accuracy with which we 
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Figure 4. Purity (top panel) and completeness (bottom panel) 
rates as a function of redshift for different dark matter halo masses 
(M^) and total stellar mass (M^^). The relationship between 
these two mass measurements is computed directly from the sim¬ 
ulation and section §4.3 is devoted to it. We find that the purity 
rates are always >60% while the completeness rates decreases as 
a function of redshift, being >80% (70%) for total stellar masses 
> 1-5 X IO^^Mq (1.2 X IO^^Mq) up to redshift 0.8. The 
selection function extracted from this analysis can be found in 
Tabled 


Table 1. Selection functions for the ALHAMBRA Survey esti¬ 
mated from the full ALHAMBRA light-cone mock catalogue 



z 

Mcl 

(Mo) 

Aql 

Mh 

(Mo) 

Level 1 (80%) 

0.2 - 0.8 

0.8-0.95 

> 0.95 

> 1.5 X 10i2 

> 1.8 X 1012 

> 51.5 

> 60.8 

> 6 X IQi^ 

> 1 X IQii 

Level 2 (70%) 

0.2 - 0.85 

0.85 - 0.95 

0.95 - 1.2 

> 1.2 X 1012 

> 1.5 X 10i2 

> 1.8 X 10i2 

> 41 

> 51.5 

> 60.8 

> 3 X lOi^ 

> 6 X IQi^ 

> 1 X IQii 

Level 3 (50%) 

0.2 - 0.8 

0.8-0.95 

> 0.95 

> 7.9 X IQii 

> 28.6 

> 1 X IQi^ 


are able to recover the redshift of the clusters from the mock 
catalogue. In Fig. we show the halo input redshift ver¬ 
sus the recovered cluster redshift for those structures that 
are recovered. We find an excellent agreement between both 
redshifts, achieving a dispersion, ctnmad = 0.0062, almost 2 
times better than the mean p hotometric redshift accuracy of 
the ALHAMBRA survey (see lMolino et al]|2014l and section 
§4.1). 


Figure 5. Purity (top panel) and completeness (bottom panel) 
rates as a function of redshift for different dark matter halo masses 
(M^) and richnesses (Acl)* Th® relationship between these two 
mass measurements is measured directly from the simulation and 
section §4.3 is devoted to it. The purity rates are always >60% 
while the completeness rates decreases as a function of redshift, 
being >80% (70%) for Acl > 51.5(41) up to redshift 0.8. The 
selection function extracted from this analysis can be found in 
Table [H 



0.2 0.4 0.6 0.8 1.0 1.2 


Figure 6. Colour-coded density map of the cluster recovered red¬ 
shift versus input halo redshift for the matched detection in one 
of the ALHAMBRA mock catalogues. The dispersion obtained 
in the cluster redshift measurement is (Tnmad ~0.0062, almost 2 
times better than the photometric redshift precision of the pho¬ 
tometric redshift of the survey. 
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log(Mj (Mg) log(Mj (Mg) 

Figure 7. Density maps of the dark matter halo mass versus 
total stellar mass for the matched haloes in the ALHAMBRA 
mock catalogue for four different redshift bins. The solid line refers 
to the linear fit of the two quantities. The vertical dotted line 
shows the mass limit to which the fit is performed. The dispersion 
measured as the standard deviation between these two variables 
is displayed, for each different redshift bin. 


4.3 Dark matter halo mass - optical richness 
calibration 

The calibration of the mass - observable relation for a clus¬ 
ter finder needs to be well understood in order to accom¬ 
plish a realistic translation of the mass. This is particularly 
cruc ial for cosmologic al purposes with galaxy cluster counts 
fe.g. lRozo et al.lf2009ll . We have measured the accuracy with 
which we can calibrate cluster and group masses with multi¬ 
ple medium-band photometry mimicking the ALHAMBRA 
data (section §4.1). 

In Fig. [3 we present the relation between the dark mat¬ 
ter halo and total stellar mass richness for the matched out¬ 
put detections in the simulations for different redshift bins. 
These relations have been used to calibrate the observable 
total cluster stellar mass in FigU) 

We now have calibrated the dispersion obtained given 
a particular richness (e.g or Aql). To do this, we have 
used a Monte Carlo approach. For each richness value, we 
have sampled randomly 10000 times all the possible halo 
mass values available and obtained a mean value and scatter. 
The results are shown in Fig. [8] for different redshift bins. 
We also quote the mean scatter, , down to the mass 

limit we reaching for each redshift bin. 

We find that the stellar mass has a very tight disper¬ 
sion in recovering masses ~ 0.27 dex) down to the 

limit of Mh ~ 3 X 10^^ Mg. This value refers to the standard 
deviation between the two variables and it is similar to the 
values that present broad-band surveys have found for other 
optical proxies such as A 200 , with the exception that we are 
able to sample 2-3 orders of magnitude lower limit in mass 



11.8 12.4 13.0 12.4 13. 

<log(M*L)>(MQ) <log(M*L)>(Mo) 


Figure 8. Average logarithm of the halo cluster mass as a func¬ 
tion of the average logarithm of the total stellar mass for different 
redshift bins. The solid line refers to the linear fit of the two quan¬ 
tities. The dotted line shows the mass limit that we are complete 
and to which the fit is performed. The dispersion measured as 
the standard deviation between these two variables is displayed, 
for each different redshift bin. This dispersion becomes compara¬ 
ble to what has been found in other optical broad-band surveys 
down to 2-3 orders of magnitude higher mass limits and to other 
non-optical mass estimators such as X-ray luminosity or Ygz- 

iIRozo et al.ll200sl : iHilbert fc Whitell20inl : lAndreon fc Bergi 
I2OI2II . In addition, this value is comparable to the precision 
that other non-optical proxies, such as X-ray luminosity or 
the total integrated SZ signal over the clus ter, Vsz, are ob- 
taining for a more extended mass range le.g iRozo et al.ll2014l 
and references herein). This result really demonstrates the 
enormous potential that multiple narrow and medium-band 
filters surveys have for calibrating cluster and group masses 
with excellent precision. 

We have performed the same analysis with the Acl pa¬ 
rameter, obtaining a very similar behaviour for the Mql, 
and a very similar dispersion. In order to show the proximity 
of these two parameters, we display in Fig. the relation 
between the Acl and the Mql parameters. As expected, 
both parameters show a very tight linear relation with very 
low scatter (0.049). Hence, we choose to work with the M^l 
hereafter. 


5 ALHAMBRA OPTICAL DETECTIONS 

We have applied the BCF to the ALHAMBRA survey with 
the specifications mentioned before. As in the simulations, 
we have pre-selected the galaxies in the catalogue with 
odds > 0.5/(l -I- z). We have also removed from the cat¬ 
alogue those galaxies with saturated flags, SaturJ'lag = 1 
and high stellar indicator indices (Stellar J'l aQ > 0.7). Be¬ 
sides, we have used the masks determined in fMolino et al.l 
(l2014h . 
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log(M*L)(M0) 


Figure 9. Density map of the logarithm of the versus the 

Acl parameter. We display the dispersion measured as the stan¬ 
dard deviation between these two variables. As expected, there is 
a very tight relation linear between these two parameters, making 
their use almost irrelevant for calibrating the cluster halo mass. 


We have set the three different thresholds determined 
in Table [T] The number of cluster and groups correspond¬ 
ing to different levels of purity and completeness are 176, 
359 and 365 for Level 1(>80%), Level 2(>70%) and Level 3 
(>50%) respectively. Note that even if Level 3 goes deeper 
in mass, it is more restrictive in redshift than Level 2 and 1, 
in order to achieve the 50% completeness and purity rates. 
We have visually checked the Level 1 and 2 detections and 
eliminated eleven of these detections (3%) that were cen¬ 
tered on a flawed part of the image, obtaining a final sample 
of 171 and 348 structures for the Level 1 and 2, respectively. 
We have also included five additional structures detected 
at 2 < 0.2 for inspection purposes. Then, for the Level 1 
detections, we find ~ 66.2 clusters and groups per square 
degree in ALH-2 to 8 fields within 0.2 < 2 0.95, while 

for the Level 2 detections, the numbers increase to 125.2 
detections per square degree within 0.2 < 2 ^ 1.2. 

We include a list of the detections in Level 1 and 2 
in Table [2] and the complementary detections to the Level 
3 in Table (3] The meaning of the columns are the follow¬ 
ing. The first column sets the name of the cluster or group 
detected, the two following columns are the cluster center 
coordinates, the fourth column is the galaxy cluster red- 
shift. The fifth and sixth columns refer to the total stellar 
mass, Mql, and the Aql parameter respectively. In Table 
(2] we have set a last column indicating the level to which 
the detection belongs (Level 1 or 2). The catalogues will be 
available in the electronic version of the journal and can be 
also found onlin^ together with a collection of colour images 
of the cluster^0. 

^ http://bascaso.net46.net/ALHAMBRA_clusters.html 
® http://alhambrasurvey.com/ 


For the purpose of this work, we will use the detection 
in Level 2, as a compromise between high purity and com¬ 
pleteness. 


5.1 Comparison with other studies 


One of the most attractive features of the ALHAMBRA sur¬ 
vey is the number of different surveys with overlapping data, 
prov iding an excellent way to deal with the cosmic variance 
fe.g. iMolino et aT]|2014l : iLopez-Saniuan et al.|[20l3i . In ad¬ 
dition, we can take advantage of early searches on these 
fields to cross-correlate them with our detections and quan¬ 
tify the degree of agreement between different datasets and 
methodologies. 

From all the surveys extending over the seven ALHAM¬ 
BRA fields, one of them has been widely studied for cluster 
purposes in the literature: the COSMOS survey, overlap¬ 
ping with the ALH-4 field. Up to date, several work has 
provided us with catalogues of large-scale structure or clus- 
ter/group detections in the COSMOS survey in the optical. 
IScoville et al.l (l2007bli detected very large-scale structure up 
to redshift 1.1 in th e COSMOS f ield u sing photometric red- 
shifts. In addition, lOlsen et al.l ll2007l i detected clusters in 
the CFHTLS-Deep Survey, also overlapping w ith the COS¬ 
MOS field. Later on, iBellagamba et aH ll201ll ') used an opti¬ 
cal and weak leasing search to dete ct galaxy clusters in the 
COSMOS field. Two recent works, ICastignani et all (l2014l i 
and IChiang et all (l2014li have provided a list of high redshift 
(2 > 1) clusters and proto-clusters detections. The former 
is based on detecting over densities around radio galaxies to 
detect high-redshift (2 1 — 2) galaxy clusters and the lat¬ 

ter, used the photometric redshift obtained from the combi¬ 
nation of the optical and IR data to detect very high-redshift 
proto-clusters (2 ~ 2 — 7). 

In addition, there are some more works that have 
been detectin g grou ps and clusters with other techniques: 
iKnobel et al.l ll2009li used a combination of a Friends-of- 
Friends (FoF) and Voronoi-Delaunay Method (VDM) al¬ 
gorithms to detect clust ers optically using the z COSMOS 
spectroscopic sample and iFinoguenov et al.l (120071 1 have per¬ 
formed an extended search of galaxies in clusters in the 36 
XMM-Newton pointings on the COSMOS field, obtaining 
calibrated weak lensing masses dLeauthaud et ai]l2010l l. 

We ha ve compared our detecti ons with th e optical de¬ 
tectio ns by IBellagamba et al.ll201ll . (Bll) and I Olsen et al.l 
I 2 OO 7 I (007) catalogues since they search for individual de¬ 
tections in the same redshift range as we do. We have also 
performed a compariso n with the X-ray detections found by 
iFinoguenov et al.l |2007l (F07) and with the spectroscopic de¬ 
tections found by Knobel et al. I I 2 OO 9 I (K09) in zCOSMOS, 
only considering those clusters with M > 3xlO^®M0.In Fig. 
[M we show the spatial distribution of all these detections 
in the ALH-4 field, consisting of 4 separated regions. For 
illustration purposes, we have included the detections found 
by F07 and K09 within 1 x IO^^Mq Si M ^ 3 x lO^^M©. 

We do not expect an exact distribution of the detec¬ 
tions, since their selection functions are built from differ¬ 
ent surveys with different depths and sets of data and us¬ 
ing different methods, with different systematics. However, 
we find a good agreement between some of the structures 
found in the field by eye. For instance, it is well-known that 
the COSMOS field has two main large-scale structures (e.g. 
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Table 2. Clusters and groups detected in the ALHAMBRA Survey with Level 1 and 2 of completeness and purity 


Name 

a(2000) 

(deg) 

<5(2000) 

(deg) 

Zest 

(IQI^Mq) 

Acl 

Level 

ALH0229.23+0108.13 

02:29:23.30 

+01:08:13.20 

0.24 

1.46 

46.62 

L2 

ALH0229.55+0108.37 

02:29:55.15 

+01:08:37.32 

0.25 

1.64 

52.46 

LI 

ALH0229.40+0108.47 

02:29:40.15 

+01:08:46.68 

0.26 

1.33 

44.27 

L2 

ALH0227.13+0110.25 

02:27:12.98 

+01:10:25.32 

0.38 

1.35 

41.73 

L2 

ALH0230.18+0105.18 

02:30:18.00 

+01:05:18.24 

0.50 

1.63 

50.26 

LI 

ALH0226.60+0107.2 

02:26:59.98 

+01:07:01.56 

0.39 

1.60 

49.31 

LI 

ALH0226.56+0103.24 

02:26:56.42 

+01:03:23.76 

0.60 

1.55 

48.60 

LI 

ALH0226.46+0108.1 

02:26:46.46 

+01:08:00.60 

0.38 

1.25 

39.17 

L2 

ALH0228.16+0110.6 

02:28:16.27 

+01:10:06.24 

1.17 

2.00 

70.75 

L2 

ALH0230.25+0103.50 

02:30:25.03 

+01:03:50.40 

0.58 

1.36 

40.86 

L2 


Table [ 2 ] is available in the online version of the article. A portion is shown for illustration. 


Table 3. Clusters and groups detected in the ALHAMBRA Survey with Level 3 of completeness and purity not included in Table 
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Name 

a(2000) 

(deg) 

<5(2000) 

(deg) 

Zest 

(1O12M0) 

Acl 

ALH0229.54+0106.3 

02:29:54.07 

+01:06:03.24 

0.26 

1.02 

33.30 

ALH0229.10+0115.13 

02:29:10.27 

+01:15:12.60 

0.23 

0.81 

24.74 

ALH0227.9+0103.46 

02:27:09.31 

+01:03:46.08 

0.36 

1.16 

36.38 

ALH0229.5+0108.1 

02:29:05.38 

+01:08:01.32 

0.38 

1.08 

34.87 

ALH0228.16+0114.46 

02:28:15.62 

+01:14:46.32 

0.39 

0.79 

24.38 

ALH0227.56+0104.22 

02:27:55.66 

+01:04:22.08 

0.39 

0.90 

28.11 

ALH0230.14+0112.47 

02:30:14.04 

+01:12:46.80 

0.69 

0.92 

31.57 

ALH0226.47+0112.49 

02:26:47.42 

+01:12:49.32 

0.64 

0.89 

28.39 

ALH0228.56+0039.10 

02:28:55.75 

+00:39:10.44 

0.26 

0.91 

27.77 

ALH0228.32+0044.37 

02:28:31.82 

+00:44:37.32 

0.22 

0.98 

30.18 


Table [ 3 ] is available in the online version of the article. A portion is shown for illustration. 


IScoville et al.l[2007al : lMolmo et aT]|20l4 l. which are basically 
recovered by all the works. We also find a basically ‘empty’ 
subfield, in agreement with all the other methods. We anal¬ 
yse in more detail the general level of agreement below. 

We have first matched the ALHAMBRA detections to 
the four different studies in order to study the level of agree¬ 
ment of the different sampl es. We have used an analogous 
FoF algorithm described in lAscaso et al. I (l2012l i. which we 
summarise here briefly. We first make a list of a friends of 
friends of every candidate, where ‘friends’ are defined as all 
the detections found within 1 Mpc. Then, we restrict this 
candidate list to those ‘friends’ whose photometric redshift 
satisfy Eq. O Finally, we select the closest ‘friend’ as the 
best match. 

We have obtained 86.36%, 84.61%, 100% and 73.30% 
of agreement with the Bll, 007, F07 and K09 samples 
respectively. Note that these values confirm the complete¬ 
ness rates found in the simulations. In addition, the accu¬ 
racy with which we recover the main redshift of the clus¬ 
ter with respect to the redshift found by the four works is 
0.021 ± 0.032 (Bll), 0.130 ± 0.167 (007), -0.007 ± 0.0095 
(F07) and 0.018 ± 0.018 (K09). 

We have also examined the agreement between the de¬ 
tections found in other work with respect to the ones found 
in this study. In Fig. [TT] we show the distribution of total 
stellar mass of the clusters detection in ALH-4. The shaded 
area histogram shows the distribution of the total stellar 


mass of those detections that have a counterpart in each of 
the different catalogues. We notice a very good agreement of 
the more massive structures (Mh > 1 x equivalent 

to Mcl > 1-8 X IO^^Mq) for K09 (>64%) and for 007 and 
F07 (>70%) while a departure of the distribution at lower 
masses. As for the comparison with Bll, we do not find 
a particular better agreement at high or low-masses being 
within the whole mass range ~50%. 

Similarly, in Fig. 1121 we show the distribution of redshift 
of the detections in ALH-4, where the dotted line and shaded 
area shows the distribution of the redshift only for those 
detections that have a counterpart in each of the different 
catalogues. For the optical works, Bll and 007, we find 
a good agreement (>70%) up to redshift 0.7 and we find 
a departure of the distribution at higher redshifts. For the 
X-ray group catalogue by F07, we find an agreement of ~ 
60% up to redshift 0.7, decreasing slightly at higher redshift. 
Finally, for the spectroscopic group sample by K09, we do 
find a very good agreement up to redshift 0.5 (>80%) and 
at redshift higher than 0.85 (>65%), whereas the agreement 
is worse within 0.5 ^ a ^ 0.85 (>43%). 

In addition to ALH -4, ALH-2 also overlaps with DEEP2 
and iGerke et al.l ll2012ll have performed a group search using 
the VDM. Their catalogue contains several groups with two 
or more members over the redshift range 0.65 < ^ < 1 . 5. We 
have compared our detections with the iGerke et al.l |2012| . 
(G12) group catalogue, restricted to detections with at least 
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Figure 10. Spatial distribution of the cluster and group detections in the ALH-4 field. The black plus symbols (and solid circles) and 
square symbols (and dotted circles) refer to the detections found in this work following the Level 1 and 2 selection respectively. The 
red diamonds and blue triangles refer to the optical detections found by Bellagamba et al. 2011 and Olsen et al. 2007 respectively. The 
green asterisk (and solid circle) and green crosses (and dotted circle) allude to the X-ray detections found by Finoguenov et al. 2007 with 
estimated masses M ^ 3 X IO^^Mq and 1 X IO^^Mq < M < 3 x IO^^Mq, respectively. Finally, the magenta asterisks (and solid circle) 
and pink crosses (and dotted circle) make reference to the spectroscopic detections found by Knobel et al. 2009 with estimated masses 
M ^ 3 X 10^® Mq and 1 X 10^® Mq < M < 3 x IO'^^Mq, respectively. The size of the circle corresponds to a 500 Kpc radius sphere at the 
redshift of the cluster. The solid lines define the four discontiguous fields of ALH-4 (see Fig. A.l in Molino et al. 2014 for a description 
of the geometry of the survey). 


three members. In Fig. 1131 we show the spatial distribution 
of all these detections in the stripe of the ALH-2 field (0.25 
deg^) which overlaps with DEEP2. Both samples have been 
restricted to the redshift range 0.65 < jz < 1 in order to 
compare complete samples in a similar range of redshift. 


We see an overall agreement between some of the struc¬ 
tures found in G12 and in our sample. Indeed, we find that 
60% of the structures that G12 finds are recovered with a 
main redshift difference of za — zq = —0.001 ± 0.021, be¬ 
ing za and zg, the redshift of this work and G12, respec- 
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Figure 11. Total stellar mass distribution for the ALH-4 field 
(solid line) overplotted with the total stellar mass distribution of 
the ALH-4 subsample that has a counterpart in each of the consid¬ 
ered works (dotted line and shaded area; Bellagamba et al. 2011 
optical/weak-lensing, Olsen et al. 2007 CFHTLS, Finoguenov et 
al. 2007 X-ray, and Knobel et al. 2009 zCOSMOS samples re¬ 
spectively). The more massive structures are well matched for 
the 007, F07 and K09 samples, finding a departure in the detec¬ 
tions at lower masses. For Bll, the rate of agreement is constant 
within the mass range. 


lively. This agreement increases to 75% if we consider those 
structures with at least five members, obtaining a redshift 
difference of za — zq = 0.001 ± 0.0015. 

We have also noticed some other structures detected in 
this work that are not recovered by G12. In order to inves¬ 
tigate which kind of detections G12 are recovering from our 
sample, we show in Fig. 1141 the total stellar mass distribu¬ 
tion of our restricted detections and overplotted, the stellar 
mass distribution only for those structures with a counter¬ 
part in the G12 sample. We see that both distributions agree 
very well finding a general agreement of 85%, being slightly 
higher (>93%) for high-mass clusters (> 7.0 x IO^^Mq). We 
then confirm that the BGF applied to ALHAMBRA-like sur¬ 
veys is able to recover low-mass groups in a similar way as 
spectroscopic surveys. 

Complementary, ALH-3 and -8 also overlap with the 
Sloan Digital Sky Survey (SDSS). Nevertheless, the limit¬ 
ing area of ALHAMBRA restricts its detection to redshift 
>0.2, whereas the bulk of the SDSS detections are below 
this redshift. 

These comparisons suggest that the optical methods are 
able to recover the redshift of the cluster similarly to the 
BGF, and with a similar precision as the spectroscopic sam¬ 
ple, at least at lower (z < 0.5) and higher redshift {z > 0.65). 
On the other hand, the X-rays and spectroscopic methods 
seem to estimate the cluster mass similarly to the BGF 
(70%), whereas two out of the three optical methods trace 
similarly the mass distribution obtained with the BGF. 



Figure 12. Redshift distribution for the ALH-4 field (solid line) 
overplotted with the redshift distribution of the ALH-4 subsample 
that has a counterpart in each of the different works considered 
(dotted line and shaded area; Bellagamba et al. 2011, Olsen et al. 
2007, Finoguenov et al. 2007 and Knobel et al. 2009 respectively). 
We find an overall good trace of the redshift range up to redshift 
0.7 for Bll, 007 and F07 finding departures at higher redshift 
range. The behavior is different for the K09 sample, finding a 
good agreement at z<0.5 and z>0.85. More details are given the 
text. 


The different cluster and gronp catalogues analysed in 
this work are built from different datasets, using a variety 
of methodologies each of them carrying their own systemat- 
ics, and with different selection functions. These facts make 
the comparison difficult, as it becomes evident in Fig. 1101 
where the different spatial distributions of each dataset be¬ 
come noticeable. Notwithstanding these facts, we still find 
a general tendency in this comparison to agree at the high- 
end of the mass distribution function {Mh > 1 x 10^“^M© at 
least) and at redshifts {z < 0.7 at least), as expected from 
the simulations. 


5.2 Properties of the ALHAMBRA groups 

As mentioned in §1, low-mass group samples within a rel¬ 
atively wide area of the sky are scarce and tend to suffer 
from observational biases. In the next subsections we review 
basic properties of the detections, we provide examples of 
such detections and we compare these results with previous 
studies. 

In this work, we do not intend to perform an exhaustive 
analysis on the properties of the galactic population of these 
structures. Rather, we mean to illustrate the variation of the 
properties of the low-mass groups with respect to the high- 
mass clusters. 
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Figure 13. Spatial distribution of the cluster and group detec¬ 
tions in one of the stripes of the ALH-2 field (0.25 deg^). The 
black plus symbols (and solid circle) and square symbols (and 
dotted circle) refer to the detections found in this work restricted 
to the redshift range 0.65 < z < 1 for the Level 1 and 2 selection 
respectively. The red asterisks (and solid circle) and the red tri¬ 
angles (and dashed circle) refer to the detections found by Gerke 
et al. 2012 at redshift <1 with a minimum number of members 
of 4, and between 2 and 4 members, respectively. The size of the 
circle corresponds to a 1 Mpc radius sphere at the redshift of the 
cluster. 


5.2.1 Stellar mass distribution 

In Figs. 1151 and [T^ we show the spatial distribution of the 
detections found in each of the fields of ALHAMBRA, pro¬ 
portionally scaled to the total stellar mass and colour-coded 
by their redshift. This distribution provides a visual way to 
study the large scale structure and filaments together with 
their cosmic variance. It becomes noticeable that the AL¬ 
HAMBRA fields have very different large-scale structures. 
Fields, like ALH-6 or ALH-7 are significantly empty, com¬ 
pared with others such as ALH-4 which displays very mas¬ 
sive structures. In fact, the well- known filaments foun d in 
COSMOS at ~ 0.4, 0.7 and 0.9 dScoville et al.l[2n07bl l are 
patent in this figure. 

It is also interesting to see that some fields, like ALH- 
5 and ALH-7 have significantly less massive clusters than 
other fields (like ALH-4 and ALH-8). Finally, we notice 
filament-like structures, i.e. clustering of groups at similar 
redshift ranges. For instance, ALH-3 clearly shows an over 
density of groups at 2 ~ 0.8, ALH-8 at 2 ~ 0.5 and 2 ~ 0.65, 
etc. 

In an effort to inspect these features in more detail, we 
have plotted in Fig.[T7]the cumulative number of clusters per 
square degree as a function of redshift for the whole sample 
(black line) and each individual field. We confirm the previ¬ 
ous analysis. The presence of substructure in the COSMOS 
field (ALH-4) is evident, finding three main sharp increases 
in the cumulative function at ~ 0.4, 0.7 and 0.9. We notice 


Figure 14. Total stellar mass distribution for the stripe of the 
ALH-2 field that overlaps with DEEP2 (solid line) overplotted 
with the total stellar mass distribution of the ALH-2 subsample 
that has a counterpart in the sample of Gerke et al. 2012. A very 
good agreement is noticeable (84.84%), particularly in the more 
massive end of the distribution. 


several changes of the slope at different redshifts for differ¬ 
ent fields. In particular, in ALH-7, we distinguish a smooth 
increase in the slope of the cum ulative function between ^ 
0.55 and 0.8, already predicted bv lArnalte-Mur et al.l ll2014h . 


5.2.2 Presence and absence of a red sequence 

The existence of a tight red sequence (RS) in galaxy 
clusters has been wide l y studied in numerous works (e.g . 
iLopez-Cruz et al]l2004 iMei et aDl2006l : lAscaso et al.ll2()oR 

and references herein). The main theories of formation of 
galaxies claim that elliptical galaxies are formed at high 
redshifts in the most massive environments and evolve pas¬ 
sively since then. One of the safest probes on the ‘univer¬ 
sal’ existence of the RS down to M > 10^‘^M© and up to 
moderate redshifts (2 < 1.6) is that clusters detected with 
methods non-dependent on the RS, such as the X-rays or 
SZ, still display a well-defined RS. However, there are a 
number of observational works that have de monstrated the 
chan g e of paradigm at high (2 > 1) redshift (iBrodwin et al.l 
I 2 OI 3 I : iMeLcLaL |2014|), and lower masses (~ several times 
lO^^M©, iFinoeuenov et ^l2007f ). Star formation increases 
in these almost unexplored ranges of mass and redshift and, 
consequently, the galactic population of these structures is 
significantly bluer and of later-type. 

We have examined the nature of the RS in the detec¬ 
tions found in the ALHAMBRA survey. The galaxies at¬ 
tributed to each cluster candidate have been selected as 
those galaxies within 1 Mpc radial distance from the clusters 
centre at the redshift of the cluster, and accomplishing the 
condition set in Eq. [2] Also, in order to make a clean selec- 
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Figure 15. Spatial distribution of the detections in ALH-2, ALH-3, ALH-4 and ALH-5 fields (top left, top right, bottom left and bottom 
right panel respectively) for the Level 2 detections. The size of the circle is directly proportional to the measured total stellar mass and 
the colour scale refers to the redshift at which the redshift is located. The solid lines define the limits of each of the fields (see Fig. A.l 
in Molino et al. 2014 for a description of the geometry of the survey). The same scale applies to all the panels. 


tion, we have considered a galaxy to belong to the cluster if 
odds> 0.5/(l + z) and StellarJFlag < 0.7. 

In Fig. 1181 we show the RS of four different detections 
belonging to different increasing mass bins (from less mas¬ 
sive in the top left panel to more massive clusters in the 
bottom right panel) at a similar redshift (~ 0.7). The red 
squares refer to those galaxies classified as early-type (spec¬ 
tral template, tb, < 5) by BPZ2.0. Looking at these exam¬ 
ples, it becomes clear the difference between the least and 
most massive structures. While the richer clusters exhibit a 
well-formed red sequence, the less massive clusters display a 
few red galaxies in place. 

In order to quantify this, we have measured the photo¬ 


metric blue fraction for each cluster at a given redshift. This 
fraction has been defined as: 


where Nb refer to the blue galaxies, considered as those 
with a bluer colour than the main expected colour of the 
red sequence at its redshift minus a typical dispersion. For 
clusters at redshift <0.7, we will consider the g — i colour, 
while for clusters >0.7, we will use the i — z colour. Nt refers 
to the total number of galaxies considered to belong to the 
cluster. 

This quantity has been measured down to a fixed abso¬ 
lute magnitude. Mi = —19.6 which corresponds to i ~ 24.5, 
the magnitude limit of ALFIAMBRA at redshift ~ 1.0, to be 
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Figure 16. Spatial distribution of the detections in ALH-6, ALH-7 and ALH-8 fields (top left, top right and bottom central panel 
respectively) for the Level 2 detections. The size of the circle is directly proportional to the measured total stellar mass and the colour 
scale refers to the redshift at which the redshift is located. The solid lines define the limits of each of the fields (see Fig. A.l in Molino 
et al. 2014 for a description of the geometry of the survey). The same scale applies to all the panels. The different distributions between 
the fields is noticeable. The ALH-4 field becomes the densest and one of the most massive, while other fields as ALH-5 or ALH-7 are 
populated with less massive clusters. 


able to account for all the galaxies independent of the red¬ 
shift dimming. Also, no odds cut has been performed here. 

We have also investigated the dependence of the pho¬ 
tometric blue fraction with the cluster/group environment, 
as a function of the redshift. To do this, we have stacked 
galaxy clusters and groups in bins of redshift and we have 
computed /s as a function of cluster total stellar mass, for 
different galaxy masses ranges, which is shown in Fig llQI 

We expect the possible contamination of the non¬ 
member galaxies to be almost negligible, due to the excellent 
accuracy of the photometric redshifts of the survey, particu¬ 
larly for the bright end of the cluster distribution. In order to 


quantify this, we have computed the blue fraction of all the 
clusters in the mock catalogue given the spectroscopic red¬ 
shift and photometric redshifts separately and we have esti¬ 
mated the contamination as the absolute difference of their 
mean values for different redshift bins. As expected, the field 
contamination source of error is more than 10 times smaller 
than the poissonian errors. Both source of errors have been 
included in the error bars in Fig llOl 

As expected, the mean /s increases as a function of 
redshift both in clusters and groups for a fixed galaxy mass. 
Also, despite the small size of the sample, we observe a ten¬ 
dency between the photometric blue fraction and the clus- 
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Figure 17. Cumulative number of clusters per square degree as 
a function of redshift for the whole sample (black solid line) and 
each separated field (ALH-2, red solid line; ALH-3, blue dotted 
line; ALH-4, green dashed line; ALH-5, dotted-dashed orange line; 
ALH-6, three dotted dashed magenta line; ALH-7, long dashed 
cyan line; ALH-8, solid yellow line). The changes of the slope in¬ 
dicate the presence of clustering at a particular redshift range, 
making evident the presence of the cosmic variance. It is notice¬ 
able the departure of the ALH-4 field with respect to the other 
ones. 


Figure 18. Colour-magnitude relation of four of the detections 
obtained in the ALHAMBRA belonging to four different bins 
from the smallest (top left) to the richest (bottom right). Only 
galaxies with odds> 0.5/(l -I- z) are shown. The red squares re¬ 
fer to those with elliptical {tf, < 5) spectral type according to 
the BPZ2.0 classification. The solid lines indicates the expected 
colour-magnitude relation at the redshift of the cluster (section 
§3). The fraction of blue galaxies appears to be higher for the 
less massive structures and few early-type galaxies are already in 
place on the red sequence in the same clusters. 


ter/group total stellar mass for a fixed galaxy mass. We 
observe, with more than 2a confidence, that the low-mass 
groups have higher /b compared to more massive clusters, 
at least up to redshift 0.45. 

Interestingly, the slope of the dependence of the blue 
fraction with the cluster richness seems also to evolve with 
redshift. At high redshift, the slope seems to be flat, whereas 
it becomes steeper at lower redshift. This result would be in 
agreement with other observational results that have sug¬ 
gested two different regimes for galaxy e volution, being t he 
environment only active at low redshift (IPeng et al.ll2O10l) . 

While this is a tentative piece of evidence, the small 
size of the sample does not allow to throw a more definitive 
conclusion. Future surveys, such as the J-PAS, will provide 
hundreds of thousands of clusters and groups with even bet¬ 
ter photometric redshift precisions, increasing the statistics 
and diminishing the error bars. As a result, it will confirm 
with high confidence these tendencies. 

The aim of this work is to describe the cluster and 
group sample detected in the ALHAMBRA survey and we 
do not attempt to explore this outcome in detail here. A 
more quantitative analysis will be carried out in a separate 
paper (Dfaz-Garcfa et al. in prep). 


6 SUMMARY AND CONCLUSIONS 

In this work, we present the first release of the optical galaxy 
clusters and groups catalogue in the ALHAMBRA survey. 


According to simulations, we are able to sample the mass 
function down to 6 x IO^^Mq, up to redshift 0.8, with 
both completeness and purity higher than 80% and down 
to ~ 3 X IO^^Mq, with detection rates >70% to the same 
redshift. At higher redshift, the mass threshold increases, 
being able to detect structures with masses > 1 x IO'^'^Mq up 
to redshift 1 in the first case, and > 6 x IO^^Mq up to redshift 
1 and > 1 X IO^'^Mq up to redshift -^1.2, in the second 
case. Additionally, we have compared our detections with 
other cluster samples obtained from a variety of datasets and 
techniques, achieving a very good agreement when matched 
to our sample, confirming the completeness of our results. 
Additionally, we have shown how the detections that are not 
found in other work are mainly located at high redshift and 
low masses, confirming the reliability of the results. 

This mass limit threshold conveys an important result, 
which is the evidence that deep multi-band medium and 
narrow-band filters allow us to sample the mass function 
with reliability down to smaller masses than deep broad¬ 
band surveys. For instance, we detected galaxy clusters 
and groups in two optical broad-ban d surveys: the DLS 
llAscaso et all [20143) and the CARS (lAscaso et al ][2 oI 3) 
surveys with the same methodology, the BCF. We were 
only able to obtain mass limits of > 1.2 x 10 ^'^Mq and 
> 4 X IO^^^Mq with completeness and purity rates >70% 
and >80% respectively, even if the data was ~ 2.5 magni¬ 
tudes deeper for the DLS and similar depth for the CARS. 

A second important consequence is the fact that the 
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Figure 19. Relationship between the total cluster stellar mass 
and the photometric blue fraction for the ALHAMBRA cluster 
and group sample stacked into different redshift bins. The differ¬ 
ent lines refer to different galaxy masses bins (1 X 10® < M* < 
1 X 10^®, dashed blue line; 1 X 10^® < M* < 1 X 10^^, dotted 
red line;l X 10^^ < M* < 1 X 10^®, solid black line). The galaxy 
member candidates have been selected within 1 Mpc distance of 
the center and performing the photometric redshift cut specified 
in Eq. [21 No odds has been performed. Blue galaxies are consid¬ 
ered those with a bluer colour than the main expected colour at 
its redshift minus 0.3 magnitudes. The error bars include simple 
Poissonian statistics and field contamination errors. 


high photometric redshift accuracy that such surveys pro¬ 
vide, allows us to obtain a reliable determination of the 
membership of the cluster. This directly translates into one 
of the highest accuracies at calibrating the redshift and mass 
of the clusters reachable with optical data up to date. In 
this work, we have been able to calibrate the total stel¬ 
lar mass - dark matter halo relation with a precision of 
0.25 — 0.35 dex down to ~ 3 x 10^®Mq, which is 
very similar to what other techniques (optical. X-rays, SZ) 
have found for, at least, two orders of magnitude higher mass 
limit. Moreover, we have measured the dispersion at recov¬ 
ering the overall redshift of the cluster or group, obtaining 
a main dispersion of ctnmad ~0.006. 

Since few optical surveys have been able to sample the 
mass range to these lower limits, the evolution of properties 
of the galaxies in those groups has not been widely explored 
due to, on one hand, the inability to obtain a complete sam¬ 
ple and, on the other, due to a possible bias of selecting ‘red 
clusters’ when using other techniques related with the CMR. 
In this study, we have preliminarily looked into the overall 
characteristics of the galactic population of these new set of 
groups detected in the ALHAMBRA survey, reaching very 
low mass limit thresholds. 

In particular, we have used the optical ALHAMBRA 
group sample to report the visible increase of the fraction of 
blue galaxies in low-mass groups. Indeed, we find a signifi¬ 
cantly lower fraction of blue galaxies in 2 ; < 0.5 ALHAM¬ 


BRA groups consistent with mo re efficient enviro nmental 
quenching in the local Universe (jPeng et al.ll201ol) and in 
agreement with other results on the blue fraction in more 
massive systems (iRaichoor fc Andreon|[201^ . |2014| ). 

These results become paramount for future applica¬ 
tions of_this_^mdof_survey. In particular, the J-PAS sur¬ 
vey ([Benitez et al.l 1^0141 ) is a survey which will be starting 
in 2015 and will image 8600 square degrees with 54 optical 
narrow-bands, down to r ~ 23.5. The expected photometric 
redshift accuracy of this survey is 0.003 and therefore, the 
expected number of groups and clusters that we can detect 
with reliability amounts to higher numbers than other sim¬ 
ilar projects aiming to go deeper with fewer, broader bands 
(Ascaso et al. in prep). As a result, this kind of data will also 
allow to confirm with high significance (>10cr) the possible 
galaxy evolutionary mechanisms happening in clusters and 
groups. 

Finally, as oth er wo rks have already cl aimed 
ijLopez-Saiiiuan et akl l2014l : lArnalte-Mur et^^ l2014l : 
iMolino et al. 2014 ). the study of the distribution of the 
cluster properties in each of the different fields of ALHAM¬ 
BRA sets evidence on the cosmic variance at the level of 
the clustering of clusters. We notice striking differences 
between different fields in terms of mass and redshift 
distribution. In future work, we will explote both the large 
scale properties of these structures and the properties of 
the galactic population of these detections 
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